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*7 TITLE OF THE INVENTION 



METHOD AND APPARATUS FOR RENDERING USING SYMPLECTIC RAY 
TRACING 

BACKGROUND OF THE INVENTION 

1 . Technical Field to which the Invention Pertains 
[0001] The present invention relates to the rendering techniques which traces the 
paths of curved light rays and renders the object.More particularly, the present 
invention relates to the apparatus and method for rendering using symplectic ray 
tracing which solves the Hamilton's canonical equation by applying symplectic 
numerical method and traces the paths of light rays. 

2 . Description of the Related Art 
[0002] It is a basic technique to visualize objects for analyzing dynamical 
systems. There are many methods of visualizing dynamical systems by using computer 
graphics. Among the many methods of visualization,ray tracing is probably one of the 
most famous methods using computer graphics. 

[0003] Conventional ray tracing employs straight lines to calculate the paths of light 
rays.The equation of a straight line is given here. 



where r is a location component on the paths of light rays,a is a starting point^v is a 
direction vector and s is a parameter. 

[0004] (Equation 1) can be rewritten as follows for the components of the vectors: 



where xi is a component of vector r.By differentiating twice,(Equation2) can be 
rewritten as a form of differential equation: 

— = 0 (Equations) 

[0005] Traditional ray tracing is thus a rendering technique in which the path of a 
light ray is calculated by solving (Equation 3). 



r = a + sv 



(Equationl) 



X i = a i + svi 



(Equation2) 
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[0006] But,there are many non-linear optical phenomena such as a quantity of air 
above the ground which is a phenomena in non-homogeneously transparent 
object.People would sometimes like to visualizing the four-dimensional black hole 
space time. 

[0007] For example,in non-homogeneously transparent object, the equation for light 
rays is given by 

dndx^ dn . 

— — — + n — = —J- (Equation4) 

ds ds ds axf 

where n is the refractive index. (Equation 4) is derived from the Fermat's principle 

in geometric optics.When n is constant,(Equation 4) reduces to (Equation 3). 

[0008] In black-hole space time, the equation for a light ray is given by 

dx^ dx^ 

r- + r' = 0 (Equations) 

ds ds ds 

where ^called the Christofifel symbol, is a function used to calculate the 
curvature of space.(Equation 5) is known as the geodesic equation. 

[0009] Traditional ray tracing can not visualize the phenomena in non-homogeneously 
transparent object or the four-dimensional black hole space time. 

SUIMMARY OF THE INVEl^ION 

[00010] An object of the present invention is to visualize the phenomena in 
non-homogeneously transparent object or the four-dimensional black hole space time. 
[00011] The present invention forms the Hamilton's canonical equation by applying 
the fast automatic differentiation techniques,and applies the symplectic Eider method 
to the formed Hamilton's canonical equation,and solves the equation by practicing the 
symplectic integration. 

[00012] The present invention fixes the location of observation and fixes view screen 
of observation,practices the symplectic ray tracing by applying the fast automatic 
differentiation techniques and S3rmplectic integration described above,and acquire the 
information of the color where the light rays cross the surface of the objects,render the 
objects according to the information of the color acquired. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[00013] Fig.l is a diagram which shows the meaning of the symplectic map. 

Fig.2A is a diagram which shows the result of the calculation obtained by the 
explicit Euler method. 

Fig.2B is a diagram which shows the result of the calculation obtained by the 
symplectic Euler method. 

Fig.SA shows the values of the Hamiltonian obtained by the explicit Euler method. 

Fig.3B shows the values of the Hamiltonian obtained by the symplectic Euler 
method. 

Fig.4 shows the construction of the apparatus for rendering of the the present 
invention. 

Fig.5 is the flowchart of the rendering process of the present invention. 
Fig.6 shows the outline of the method for rendering of the present invention. 

Fig. 7 shows the example of the rendering obtained by the present invention. 
Fig.8 shows the example of the rendering obtained by the present invention. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[00014] First,the principle of the present invention is described. As descried 
above,past studies have adopted equations such as the geodesic equation,the 
Helmholtz equation,the Eilconal equation,and so on.However»from a mathematical 
point of view,all of those equations can be represented in terms of the Hamilton's 
canonical equation: 

dp. __dH 
ds dq- 

da. dH 

- - — (Equations) 
ds op. 

where qi indicates the location component, pi indicates the momentimi of qi and H is a 
Hamiltonian composed of qi and pi . 

[00015] For example,the Hamiltonian that lead to the (Equation 3) is given by 
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H = (px2 + Py2 + Pz2 )/2 

[00016] When the Hamiltonian is partially differentiated,and form the Hamilton's 
canonical equation,the equation of the straight line is obtained.In non-homogeneously 
transparent object such as a quantity of air above the ground,the Hamiltonian is given 
as 

H = n (x,y,z)(px2 + py2 + pz2 )/ 2 
[00017] In consideration of the theory of relativity,the Hamiltonian for calculating 
the paths of light rays in spherically S3m[imetrical black-hole space-time is given as 



H = 



r 



2 I r J 2 2r^ 2r^sin^e 
where (f,r,0, 0) are the components of a set of four-dimensional polar coordinates, 
(Pt^ Pr9 Pe^ Pip) denotes the corresponding momentum components^and rg is the mass 
corresponds to the mass of the black hole called the black-hole radius. 
In general, provided the Einstein equation is gM , a Hamiltonian for the light rays is 
given as 

H = ^±g'p,Pj 

[00018] As described above, the Heuniltonian of the motion of the light rays can 
usually be obtained.Therefore, various objects can be simulated simply by applying the 
Hamilton's canonical equation without implementing different equations. 
[00019] Here the S3niiplectic numerical method of the present invention will be 
described.The characteristic of the Hamilton's canonical equation is that the equation 
preserves the property of the S3rmplecticness.The symplecticness is the property of the 
symplectic map defined as follows. 

[00020] (Definition) Let(p(to),q(to)) be a point at (t=to) and (p(to+ A t ),q(to+ At)) be 

a point at (t=to+ At).For all real numbers to, the mapping (p(to),q(to))— > (p(to+ A 

t ),q(to+ A t)) is a symplectic map if and only if dq(to) A dp(to)=dq(to + A t) A dp(to+ A 

t),where dq and dq are one-forms of v£udables and A is the exterior product operator. 
[00021] Fig.l is a diagram which shows the meaning of the sjrmplectic map.The 

symplecticness means that the area (S) of the parallelogram constructed of the vector 
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dp and dp is preserved. 

[00022] The present invention applies the Hamilton's canoniced equation to any 
objects and solves the equation by applying the symplectic symplectic numerical 
method which preserves the property of- the symplecticness and traces the path of light 
ray. This method is called the symplectic ray tracing.In the S3anplectc nixmerical 
method, the result of the calculation adopting original value satisfies the definition 
above. 

[00023] The present invention is able to reduce the error occurred in tracing the path 
of the light ray by applying the symplectic integration. Therefore the phenomena in 
non-homogeneously transparent object or the four-dimensional blacl^ hole space time 
can be visualized. 

[00024] The present invention uses the symplectic Euler method which is the 
implicit Euler method as the symplectic numerical method.The symplectic Euler 
method applies the scheme shown below to the Hamilton's canonical equation. 

^<lk^i 

[00025] The performance of the calculation of the symplectic Euler method differ 
much from the performance of the esplicit Euler method which uses the scheme shown 
bellow. 

Pm = Pk 5 

[00026] Fig.2A and Fig.2B shows a very simple example comparing S3ncnplectic Eider 
method and explicit Euler method.Each curve in the (p-q) plane indicates the orbit of a 
light ray where H = pq(p+q-l). 

[00027] Fig.2A shows the result obtained by using the explicit Euler method while 
Fig.2B shows that obtained by using the symplectic Euler method.Since analytical 
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solutions of the equation for this system give closed curves, Fig.2B should give a 
correct result. 

[00028] Fig.3A and Fig.3B shows the values of the Hamiltonian that have to be 
preserved.Fig.3A shows the values of the Hamiltonian obtained by the exphcit Euier 
method.The Hamiltonian decreases as the calculation proceeds. 
[00029] Fig.SB shows the values of the Hamiltonian obtained by the symplectic 
Euler method.The values of the Hamiltoniem are preserved as long as the vibrations 
are negligible. 

[00030] In this test case,the form of the Hamiltonian is so simple that it appears 
there are no problem with the symplectic nimierical method. Ho we ver,it is hard to 
apply the Hamilton's canonical equation to a complex Hamiltonian^because forming 
the Hamilton's canonical equation by applying the sjmiplectic Euler method requires 
all the derivatives of the Hamiltonian. 

[00031] For example ,the Hamiltonian in axisymmetric space-time is given as 



Pt A 1 



V ^ 

2 2 

s 



2 2 

(Equation?) 



r„r ar.r 



p^Asin^e 2 "^p^A^'^*' 
where a corresponds to the angular momentiom of the black hole,and 

For the above Hamiltonian,Mathematica,a software for culclation, outputs the 
derivative of 0 as follows. 
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2 Pa P,rr^ cos 6 sin Oa 



p^^(r^ -r^r + a^ cos^0) cot 
(a^ + - rr J(r^ + a^ cos^ 0)^ 

p^^ cot0a^ 
(a' + - ^^•^K^'' + cos^ 0) 

cosOsinOa^ 
(r^ + a^cos'0)' 
Pr^ia^ - rr^ ) cos 0 sin 
(r' +a'cos^(0))^ 
Pj^r(a^ cos0sin0a^ 
(a^ + - rr^Xr^ + cos^ 0)^ 

p.^(r^ - rr+ a^ cos^ 0) cot 0 csc^ 0 
+ — I £ 

(a^ + cos^ 0) 

[00032] Implementation of such a complicated expression would lead to mistakes 
and a slow computation speed. 

[00033] The present invention adopts the fast automatic differentiation techniques 
and calculates the differentiation of a given Hamiltonian H automatically and forms 
the Hamilton's canonical equation. 
( Fast automatic differentiation techniques ) 

[00034] In this fast automatic differentiation techniques, a vector is defined 
(fip,q),dpf,dqf) 

For example.addition and mtdtiplication are defined as follows: 

if,d^f,d^f)xig,d^g,d^g) 

= ( fg,d^fx g+fx d^g,d^ fxg+fx d^g) 

[00035] As shown above,by defimng the overload of the operator.the operations and 
partied differentiations are calculated simultaneously.For example.by using the rule of 

the calculation shown above, the result of multiplication as follows can be obtained, 
(x2,2x) X (x,l)=(x3,3x2),while x2 times x is x3 and the differential is 3x2.That is to 
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say,when the multiplication is calculated, the differential can be obtained 
automaticsdly. 

[00036] When the value of the Hamiltonian is calculated by applying the sjmiplectic 
ray tracing, all the derivatives of the HamUtonian can be obtained 
simultaneously.Therefore,the Hamilton's canonical equation can be formed. 
[00037] For a given Hamiltonian H = pq, fast automatic difiTerentiation calculates 

(p,l,0)x(q,0,l) 

= iP^<l^d ^pxq-\- pxd ^q,d^pxq+ pxd^q) 
= (pq,lxq+ pxO,Oxq+ pxl) 



As shown above,the operations and partial differentiations are calculated 
simultaneously without truncation errors. That leads to the formation of the 
Hamilton's canonical equation. 

[00038] Fig.4 shows the apparatus for rendering which the present invention 
provides. 1 is the apparatus for rendering.il is the section of fixing the location of 
observation and view screen of observation. 12 is the section of symplectic ray tracing 
which practices symplectic ray tracing. 13 is the section of acquiring the information of 
the color which acquires the information of the color of the cross section of the light 
ray and the surface of the objects. 14 is the section of rendering which renders the 
objects based on the information of the color acquired. 15 is the section of storing the 
information of the object which stores the information of the object and of the 
Hamiltonian H which corresponds to the object, 

[00039] In the the section of symplectic ray tracing 12,121 is the section of forming 
the Hamilton's canonical equation which forms the Hamilton's canonical equation, 122 
is the section of practicing symplectic integration which practices symplectic 
integration by applying the symplectic Euler method to the Hamilton's canonical 
equation. 

[00040] Following,the rendering process of the present invention will be described by 
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using the Fig.5 and Fig.6.Fig.5 is the flowchart of the rendering process of the present 
invention.In Fig.6 which shows the the outUne of the method for rendering of the 
present invention,the string CRL 30 in the lower left part represents the object. 
[00041] First, the section of fixing the location of observation and fixing view screen 
of observation fixes the location of observation(step Sl).For example the section fixes 
the location of the eye of the observer 31 shown in Fig.6. Next, the section of fixing the 
location of observation and fixing view screen of observation fixes view screen of 
observation 32(step S2)Jn Fig.6,view screen of observation 32 is shown in the 
middle.Next,the section of practicing the S3miplectic ray tracing 12 practices the step3 
and step4 shown below to trace the path of light ray until the light ray cross to the 
object(simplectic ray tracing). 

[00042] The section of forming the Hamilton's canonical equation 121 forms the 
Hamilton's canonical equation(step S3).That is to say,based on the Hamiltonian H of 
the object which is deribed fi*om the section of storing the information of the object 15, 
the section of forming the Hamilton's canonical equation 121 forms the Hamilton's 
canonical equation applying the fast automatic differentiation techniques. For 
example,as described above,when H=pq is calculated by applying the Fast automatic 
differentiation techniques, all the derivatives of the Hamiltonian shown below can be 
obtained and the Hamilton's canonical equation can be formed. 

(Al,0)x(q,0,l) 

= ipxq.d^pxq+ pxd^q^d^pxq-h pxd^q) 
= (pq,lx q+ pxO,Ox q+ pxl) 



[00043] And the section of practicing symplectic integration 122 practices symplectic 
. integration(step S4). That is to say, the section of practicing symplectic integration 122 

practices symplectic integration by appl3mig the sjmiplectic Buler method to the 
formed Hamilton's canonical equation. 

[00044] By practicing the symplectic integration by applying the symplectic Euler 
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method, the value of p and q can be obtained one after another and the paths of the 
light ray can be traced. 

[00045] Next,the section of acquiring the information of the color 13 acquires the 
information of the color(step S5).In particular^the section of acquiring the information 
of the color 13 acquires the information of the color of the cross section of the light rays 
and the surface of the object.And the section of rendering the objects 14 renders the 
object(step S6).And the process described above finishes.For example,when the Ught 
rays travel along the dotted curves,the string CRL is rendered upside down.33 is the 
illustration rendered on view screen of the object 32,and 34 is the illustration perceived 
by the observer's eye. 

[00046] The exsunple of the rendering of the present invention is shown below.In the 
present invention^blacl^-hole space-time is refered as a four-dimensional space with a 
black hole at the origin. For example, the Hamiltonian as follows generates the result 
shown in Fig.7. 



r 



\ 2 /" rA_2 _2 



P, 



1-^ 



V 



r j 2 2r^ 2r^sin^e 
This phenomenon is Imown as the gravitational convex lens effect. 
[00047] The Hamiltonian of (Equation?) which is in an sodsymmetric space-time 
generates the result shown in Fig.8. 
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